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The Jahn-Teller effect is a fascinating and ubiquitous phenomenon in modern quantum physics 
and chemistry. We propose a class of oxides with melilite structure Ba 2 TGe 207 (T=V,Ni) where 
Jahn-Teller distortions are the main responsibles for the appearance of electric polarization. At the 
hearth of the proposed mechanism lies the lack of inversion symmetry displayed by tetrahedrally 
coordinated transition-metal ions, which allows for the condensation of polar Jahn-Teller distortions, 
at odds with octahedral coordination typical of conventional ferroelectric oxides with perovskite 
structure. Since the noncentrosymmetric local environment of transition-metal ions also activates the 
proposed spin-dependent hybridization mechanism for magnetically-induced electric polarization, 
proper multiferroic phases with intrinsic magnetoelectric interaction could be realized in this class 
of low-symmetry materials. 

PACS numbers: 77.80.-e, 71.70.Ej, 75.85.-l-t 


I. INTRODUCTION 

Jahn-Teller (JT) distortions represent a universal 
mechanism by which spontaneous symmetry breaking 
may occur in condensed-matter systems. Colossal mag¬ 
netoresistance in manganites has been explained invok¬ 
ing an essential role of the Jahn-Teller effecti, which 
has been also invoked in several high-temperature oxide 
and fullerene superconductors^i^. Even though JT dis- 
torsions are nonpolar in the perovskite ABO 3 structure 
displayed by many ferroelectric oxides, the ferroelectric 
transition in these systems has been also explained in 
terms of a pseudo Jahn-Teller (PJT) instability which 
can take place whenever the vibronic coupling between 
the ground and excited states is sufficiently strong^. The 
cooperative PJT effect in ferroelectrics has been known 
for a long time, nevertheless its possible realization in 
multiferroic perovskite oxides, showing the simultaneous 
presence of magnetic and ferroelectric ordering, has been 
only recently suggested^^— .Interestingly, such mechanism 
defies the empirical exclusion rule, according to which 
proper ferroelectricity and magnetism should be chemi¬ 
cally incompatible and mutually exclusive^. 

On the other hand, the quest for sizeable magne¬ 
toelectric (ME) effects has been pursued in the last 
decade mainly among magnetically induced improper 
ferroelectrics, where the microscopic ME interaction is 
expected to be intrinsically large. In this respect, a 
very interesting class of materials has been recently ob¬ 
ject of intense research activity, suggesting a local ori¬ 
gin of ME coupling in low-symmetry crystals through 
a spin-dependent hybridization mechanisnti^rJ^. Es¬ 
sentially, the lack of inversion symmetry in the point 
group of the transition-metal ions may drive the ap¬ 
pearance of local dipoles via an anisotropic hybridiza¬ 
tion, modulated by the atomic spin-orbit coupling, be¬ 
tween the transition-metal ions and the surrounding oxy¬ 
gens. Materials belonging to the melilite family, with 
general formula A 2 TB 20 j alkaline metal, tran¬ 


sition metal, 5+^ =Si, Ge), consisting of B 2 O 7 dimers 
linked by TO 4 tetrahedra, represent an interesting play¬ 
ground, where the local properties of cations with 
tetrahedral coordination could be explored in details. 
Among these, several compounds have been already syn¬ 
thesized, i.e. Ba 2 rGe 207 (T =Mn,Co,Cu) ^^d^ and 
A2CoSi207(A =Sr,Ca)iiii^. Due to the lack of inversion 
symmetry and to the quasi-two-dimensional character of 
their magnetic interactions, these materials have been 
predicted to host peculiar incommensurate magnetic spi¬ 
ral ordering!^ and skyrmion excitations^^, multiferroicity 
and highly non-linear magnetoelectric respons o^^'^^’^^'^^ , 
as well as magnetochiral^^ and a giant directional dichro- 
ism in resonance with both electrically and magnetically 
active spin excitations^^. The latter phenomena have 
been mostly explained in terms of the local ME interac¬ 
tion discussed above; specifically, the multiferroic phase 
observed in Ba 2 CoGe 207 has been shown to have an im¬ 
proper origin, the ferroelectric polarization being induced 
by the magnetic order. 

In this respect, a yet unexplored possibility is to con¬ 
sider JT instabilities and possible proper ferroelectric 
transition for T ions showing degenerate electronic states 
in tetrahedral symmetry. The lack of inversion symmetry 
of the local tetrahedral environment allows for odd ionic 
displacements to couple with the degenerate electronic 
states and, therefore, for symmetry-allowed appearance 
of local dipole moments. In principle, this could imply 
proper multiferroicity with intrinsically large ME cou¬ 
pling in the family of melilite oxides. We notice that, 
despite the possible occurrence of polar JT distortions 
in crystals without local center of inversion has been 
known for a long time and experimentally revealed, e.g., 
in rare-earth oxides (such as DyV 04 ) or in praseodymium 
compouds (such as PrGla)^^— , so far its relevance in 
the field of multiferroics has been overlooked. Based on 
these premises, we first discuss in Section m the Jahn- 
Teller problem relevant for the tetrahedral units TO 4 
appearing in the melilite crystal. On the basis of gen- 
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eral symmetry considerations, we show that polar JT 
distortions may locally develop in this structure when 
JT-active transition-metal ions are considered. Density 
functional theory (DFT) calculations are then used to 
propose novel melilite oxides where such polar JT effects 
can appear, possibly leading to ferroelectric, antiferro- 
electric or ferrielectric phases, as discussed in Sections 

Unnvl 

II. POLAR JAHN-TELLER EEEECT AT 
TETRAHEDRAL SITES OF MELILITE 
STRUCTURE 


Due to the quasi-layered structure of the melilite 
structure, each r04 appears to be slightly compressed 
along the c axis, thus belonging to the tetragonal D2d 
group. As a consequence of the tetragonal crystal 
field, transition-metal d states are split into nondegener¬ 
ate dz2,dx2-y2, two-fold dyz,dzx and nondegenerate dxy 
states. The corresponding Jahn-Teller problem for de¬ 
generate dyz , dzx is then expressed as A 0 (oi -I- 02 -I- -b 
62)- By introducing the standard symmetrized displace¬ 
ments Qo{Ai),Qa{A2),Qi{Bi) and Q2(i?2), describing 
totally and non-totally symmetric displacements shown 
in Fig. [Ue)-(h), the vibronic matrix defined on the de¬ 
generate electronic states reads^: 


where 



( 1 ) 


Ai — FiQi + L2QaQ2 + GiQiQq 

^2 = F2Q2 + G2Q0Q2 + LiQiQa, ( 2 ) 

including both linear Fi,F2 and quadratic Gi = G[Ai x 
Bi),Li = G{A2 X Bi) couplings. Notice that the non- 
totally symmetric Qi and Q2 modes describe a nonpolar 
and polar displacement, respectively, as highlighted in 
Figs. [Ug)-(h). Upon inclusion of the elastic term iJo = 
(1/2) ATiQf, the leading distortions as a function of 

all vibronic parameters can be obtained by minimizing 
the corresponding energy £ = \ Y^iKiQ^ ± + 

Analytical solutions are found only when the quadratic 
Li couplings are neglected; however, all modes Qi of the 
full JT problem can be obtained as a function of the 
stationary electronic states, namely of the mixing angle 

through which the ground state is expressed as \tf}) = 
cos (p\yz) + sin ip\zx)^. 

When Li = 0 , two kinds of energy minima are found, 
where only one of the non-totally symmetric modes with 
Bi and B2 symmetry is activated, either Qi = 0, <52 7^ 0 
or Qi 7b 0, Q2 = 0, depending on the ratio between lin¬ 
ear vibronic and elastic energies K2F1 ^ K1F2 as shown 
in Fig. [TJa), (c) and later on labeled as case a and b, 
respectively^. In terms of the electronic wavefunctions, 
the two energy minima correspond to ipa = ±7r/4 (case 
a) and (/?{, = 0,7r (case 6), i.e. to a symmetric mixing of 
dyz, dzx states or to split levels with dyz or dzx unique 


character, respectively. On the other hand, when Li ^ 0 
an effective coupling between nonpolar and polar modes 
appears through a term {F1L2 + F2Li)QaQiQ2 in the 
adiabatic potential energy £. As a consequence, all dis¬ 
tortion modes are found to be nonzero at energy min¬ 
ima, thus implying that a polar displacement is always 
symmetry-allowed, as shown in Fig. O; at the same 
time, JT split states display an asymmetric mixing of 
dyz and dzx- Furthermore, the switching of Q2 (hence of 
the local dipole) implies a change of sign of the product 
QaQi, consistently with the trilinear term coupling the 
two nonpolar and the polar modes. As a consequence, 
the energy barrier related to the dipole switching can in 
principle be much reduced with respect to the case of a 
single-mode instability [see, e.g., Fig. [I](d)]. 




FIG. 1 . (a)-(d): Typical evolution of JT modes — shown in 
panels (e)-(h) — and energy as a function of the electronic 
angle tp for linear and quadratic couplings. Parameters have 
been chosen as K\ = K2 = A'a/2 = Kol2 = 1 eVjk?, Fi = 
0 . 3 ( 0 . 1 ) eV/A , F2 = 0 . 1 ( 0 . 3 ) eV/k in top (bottom) panels, 
corresponding to cases a and b discussed in the main text. 
Quadratic coupling constants are set to zero in left panels 
(a),(c), while Li = L2 = 0.9 eU/A^ and G\ = G2 = 0.5 
eV/k^ in right panels (b), (d). Arrows and vertical lines 
highlight the energy minima, (c)-(f): Schematic represen¬ 
tation of the symmetrized displacements in each tetrahedral 
unit, corresponding to totally symmetric displacements with 
symmetry Ai (e) and A2 (f), and to non-totally symmetric 
displacements with symmetry Bi (g) and B2 (h). 
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III. JAHN-TELLER EFFECT IN BaaNiGeaOr 


On the basis of the previous general analysis, we con¬ 
sider a hypotetical melilite oxide hosting JT-active T 
ions, where the polar JT effect previously discussed is 
expected to appear. We performed DFT calculations 
resorting to the VASP code^ within the framework of 
generalized gradient approximation GGA-PBE, while the 
electronic correlation has been also taken into account 
by using the GGA-l-17 potential^ (with 17=2, 4 or 6 
eV for T ions). Electric polarization has been evalu¬ 
ated in the framework of the Berry-phase approach^. 
We started from the experimental crystal structure of 
Ba 2 GoGe 207 belonging to P42 ito space group^ and 
fully optimized it by substituting Go with the JT-active 
ion while imposing different magnetic configura¬ 

tions of the transition-metal magnetic moments. We also 
considered V^+, whose electronic configuration is phys¬ 
ically equivalent to Ni^’*', but with occupied majority 
spins only. Indeed, results for both systems show qualita¬ 
tively the same behavior, therefore we discuss here only 
the case of Ba 2 NiGe 207 (BNGO). 

In Fig. [UJd), we show the density of states of BNGO 
decomposed into d orbital states, as obtained when a C- 
type antiferromagnetic configuration is imposed (similar 
results are obtained with other higher-energy magnetic 
configurations). When keeping the same space group 
P42 ito of the parent compound, BNGO clearly shows 
a metallic behavior, arising from the twofold degenera¬ 
cies of {dyz, dzx) orbital states which are half-filled in the 
minority spin manifold of Ni^+((i®), being 
The degeneracy can then be lifted by a JT structural 
distortion; indeed, we found that lowering the symme¬ 
try to c-centered Cmm2 space group leads to the open¬ 
ing of an energy gap in the BNGO compound, being 
Eg = 0.3 (2.3) eV at the bare GGA {GGA+U, 17=4 eV) 
level, with a corresponding total energy gain of 78 (753) 
meV per formula unit (f.u.). We notice that the main 
effect of a finite U in the GGA-|-t7 approach is to induce 
a larger gap and a larger energy gain associated to the 
structural transition. Both the Ni-0 bond lengths and O- 
Ni-0 bond angles are differentiated, as shown in Fig. [3] 
(b). Two short s (long 1) bonds can be identified, linking 
the Ni ion with upper (lower) lying oxygens in the tetra¬ 
hedral cage, being s=1.95 A and /=2.00 A (s=1.97 A, 
1=1.98 A for 17=4 eV). Furthermore, the 0*°P-Ni-0''‘™ 
bond angles /3, comprising upper- and lower-lying oxy¬ 
gens, differentiate in /3'=106.1° (107.6° for 17= 4 eV) and 
/3"=103.0° (102.5° for U= 4 eV), suggesting that the JT- 
induced distortions do not consist only in a Ni offcenter- 
ing. Indeed, the ionic displacement from the P42im to 
the Cmm2 structure can be decomposed^^ into a totally 
symmetric Fi mode, with total amplitude Qri = 0.02 A, 
and a distortional (non-totally symmetric) mode, with 
total amplitude Qts = 0.13 A. The latter displacement 
mode can be further decomposed in two modes belonging 
to Bi and B 2 symmetry representations, shown in Fig. 
me) and corresponding to those shown in Fig. [ljg),(h). 




Energy (eV) 


FIG. 2. (a) Melilite crystal structure of Ba 2 TGe 207 with 

P42im symmetry, showing layers of Ge 207 dimers linked by 
TOi tetrahedra intercalated by Ba ions. The r04 tetrahedra 
are compressed along the c axis, resulting in different O-T-0 
angles a, j3. (b) Top view of the crystal structure, highlighting 
two inequivalent r 04 tetrahedra which are rotated about the 
c axis of ±K, respectively, (c) Level structure of T'^'^ metal 
ion in the melilite crystal. In the tetrahedral environment, d 
orbitals split into lower Cg and higher t 2 g manifolds; a tetrag¬ 
onal compression further split the eg,t 2 g levels, leaving only 
twofold degenerated dyz,dzx states. If these levels are half- 
occupied, as in the case of Ni(d®) or V(d®), a JT distortion 
can take place, removing the only left degeneracy, (d) Den¬ 
sity of states (States/eV) decomposed into d orbital states 
and calculated within the bare GGA approach. “Ni*” refers 
to the parent-compound structure, belonging to space group 
P42im, while the optimized distorted structure with Cmm2 
symmetry is labeled by “Ni”. 


in perfect agreement with our previous qualitative anal¬ 
ysis. These modes display significant ionic displacements 
of O ions around the Ni cation, being (5 bj=-0.I0 A and 
(5 b2=0.05 a. Upon structural distortion, the degenerate 
yz,zx states are split into occupied O.Gljj/z) — 0.79|zx) 
and unoccupied 0.79|?/z) -I- 0.61 |zx) states, implying a 
mixing angle ip = 37.7°. From the point of view of struc¬ 
tural distortions, the Bi mode, differentiating both the 
bond angle and the bond length, appears as the largest 
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FIG. 3. Charge density plot (obtained within the GGA+17 
approach, with [7=4 eV) corresponding to the highest oc¬ 
cupied state of Ni-d state hybridizing with tetrahedral 0-p 
states, (a) in perspective view in a local xyz frame, (b) pro¬ 
jected onto the xy plane, (c) Distortion modes of tetrahedral 
O ions surrounding Ni ion, projected onto the xy plane, as 
obtained using the ISODISTORT program^S. 


JT mode; however, the electronic mixing angle appears 
closer to the value <pa=T’'/‘i of case a, suggesting that 
indeed the leading JT distortion arises from mode i? 2 , 
which is the main responsible for the d-orbital state split¬ 
ting. On the other hand, the simultaneous activation of 
both Bi and -B 2 modes, alongside the nonzero ampli¬ 
tude of the totally symmetric displacement, points to a 
non-negligible quadratic couplings Li, which is further 
confirmed by the asymmetric mixing of the dyz and 
orbital states. 


IV. MULTIFERROICITY IN JAHN-TELLER 
MELILITES 

If on the one hand the appearance of local dipoles in 
TO 4 cage can be understood in terms of JT effect, on the 
other hand the onset of an ordered ferro- (FE) or antifer- 
roelectric (AFE) phase is related to non-local interactions 
between the different tetrahedral units. Since there are 
two TO 4 tetrahedra in the parent unit cell, two possible 
structural configurations can be considered which real¬ 
ize the previously discussed local JT effect ; i) both T 
ions offcenter downward or upward (FE, Pc ^ 0, space 
group Cmm2) , ii) T ions offcenter in opposite directions 
(AFE, Pc = 0, space group P2i2i2i). As listed in Tab. 
m the AFE structure in both compounds is energetically 
more stable by 20-30 meV/f.u.. The antiferro character 
of the cooperative JT interactions can be qualitatively 
understood in terms of the local coupling between on-site 
JT distortions. In fact, the sign of the polar distortion 
Q2 is determined by the product QaQi, where the first 


nonpolar distortional mode is associated with a global 
rotation of the tetrahedral unit and the second to the 
antiphase rotation of and around the tetrahe¬ 
dral 0 axis, as shown in Figs UK f) and (g), respectively. 
The sign of Qa is opposite in the two TO 4 units, which 
are rotated by an angle ±k in the melilite crystal; on the 
other hand, Qi is expected to display the same sign, in 
order to minimize the energetic cost associated with dis¬ 
tortions of the Ge 2 07 dimers linking the JT-active units 
by inducing less asymmetric changes of the O-Ge-0 bond 
angles. The opposite sign of QaQi thus may explain the 
observed antiferro configurations of local dipoles associ¬ 
ated with the Q 2 local modes. However, one can resort 
to a different JT ion occupying T sites, such as V^+. In¬ 
deed, when combining V and Ni ions, each carrying a 
different electric dipole, a ferrielectric configuration can 
be expected, as we show in Tab. H] On the other hand, 
the substantially local mechanism leading to the forma¬ 
tion of electric dipoles suggests that ferroelectric behavior 
could be in principle attained when doping melilite ox¬ 
ides with JT-active ions; indeed, good-quality crystals of 
Ba 2 Gui_a;Ni 2 ,Ge 207 have been successfully synthesized 
up to Ni concentrations of x = 0.5, their structural and 
magnetic characterization being under wa y^^’^^ . 



AE (meV/f.u.) 

Pc {iiC/cra^) 

BNGO 

PE 

775.9 

- 


FE 

22.9 

1.18 


AFE 

0 

0 

BVGO 

PE 

675.3 

- 


FE 

30.3 

1.90 


AFE 

0 

0 

B(V,N)GO 

PE 

339.7 

- 


strong FE 

189.7 

-10.85 


weak FE 

0 

0.70 


TABLE I. Relative GGA-|-t/ (t/=4 eV) total energies and 
bulk polarization for paraelectric metallic (PE), FE and AFE 
configurations in BNGO, Ba 2 VGe 207 (BVGO) and an or¬ 
dered half-doped Ba 2 Vo. 6 Nio, 5 Ge 207 (B(V,N)GO) hypoteti- 
cal compound, where the parallel(antiparallel) configuration 
of inequivalent electric dipoles results in a strong(weak) FE 
phase. 

As for the magnetic properties, we found that both 
Ni and V oxides display an antiferromagnetic ground 
state, and therefore they can be considered as proper 
multiferroic materials. The in-plane exchange constant 
Jab appears to be the dominant antiferromagnetic in¬ 
teraction, while the out-of-plane exchange Jc changes 
from ferro- to antiferromagnetic interaction when Ni is 
replaced by V, leading to G-type and G-type AFM con¬ 
figurations, respectively. In both cases, the estimated ex¬ 
change anisotropy is rather strong, being Jc/ Jab ^ 0.05, 
thus displaying the quasi-2D character observed in known 
parent compounds. On the other hand, the most impor¬ 
tant contribution to the local spin-dependent hybridiza¬ 
tion mechanism mediating the reported large ME inter¬ 
action in melilite oxides has been shown to arise from 
dyz , dzx orbital states of the T ionsi^. Since these states 
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are partially filled in JT melilites, a significant spin- 
dependent modulation of charge density at lower- and 
upper-lying O ions through the asymmetric pd hybridiza¬ 
tion is expected to further contribute to the local elec¬ 
tronic polarization, i.e., to mediate a ME interaction. 
Due to the partial occupancy of the active d states of 
T ions, such ME interaction can be in principle of the 
same order of magnitude of that predicted for Co and 
Cu compounds, if not larger. Even though the relatively 
weak exchange interactions are responsible for very low 
Neel temperatures, signatures of such ME coupling are 
expected to be visible even in the paramagnetic phase, 
due to the localized nature of its microscopic origin. 

V. CONCLUSIONS 

By combining symmetry-based JT analysis and DET 
calculations, we put forward a novel mechanism for mul- 
tiferroicity and theoretically predict the possibility of re¬ 
alizing proper EE and AFE phases in melilite oxides. We 
have shown that JT instabilities of noncentrosymmetric 
units can lead to polar distortions and possibly mediate 
a (anti)ferroelectric transition of proper character which 
could be signalled by large anomalies in the dielectric sus¬ 
ceptibility. The cooperative JT origin of such structural 
transition is also known to lead to enhanced electrostric- 


tive and magnetostrictive responses, that are expected to 
develop even if the predicted JT distortions display a dy¬ 
namical, rather than static, character—. Remarkably, the 
JT polar distortions can be realized in the presence of a 
magnetic phase, thus circumventing the empirical exclu¬ 
sion rule between magnetism and ferroelectricity usually 
invoked for perovskite oxides. On the other hand, due 
to the proper nature of the predicted structural distor¬ 
tions, they are expected to take place at higher tempera¬ 
tures than those typically found for magnetically-induced 
improper ferroelectrics, at the same time displaying an 
intrinsic ME interaction of local origin that could be de¬ 
tectable even above the magnetic transition temperature. 
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